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^ ; Abstract 

One-nucleon emission electron scattering experiments are studied with a model that considers 
' short-range correlations up to the first order in the number of correlation lines. The proper normal- 

. ization of the many-body wave functions requires the evaluation of two- and three-point diagrams, 

the last ones usually neglected in the literature. When all these diagrams are included the effects 
of the short-range correlations are rather small. The results of our calculations are compared with 
^ ' experimental data taken on ^^O. 

c ■ 



X 



PACS number(s): 21.10.Ft, 21.60.-n 



1 Introduction 



One-nucleon knock out experiments induced by electromagnetic probes have been, and still are, an 
important tool to study the structure of the atomic nucleus ||l|, |2| . The main use of these experiments 
has been the investigation of the ground state single particle structure of nuclei Q. For this reason, 
they have been done, in the great majority of cases, in the quasi-elastic region, in order to minimize the 
effects of the collective nuclear excitations. Nucleon emission in the quasi-elastic region is dominated 
by the direct knock out of the detected nucleon while the other nucleons behave like spectators. 
Many-body effects can be treated as corrections to this basic process. 

These (e,e'p) experiments have been analyzed in terms of mean field models with great success. 
The models describe the target nucleus ground state by using a purely real mean field, and the nuclear 
excited states in terms of one-particle one-hole excitation with the particle moving in a complex optical 
potential. The optical potential is usually taken from a fit to proton elastic scattering data off the A-1 
nucleus. This approach reproduces rather well the behavior of the cross section in terms of the so-called 
missing momentum. What remains, however, is the problem of the overestimation of its magnitude. 
The quenching factor necessary to reproduce the data has been related to the spectroscopic factor 
which nowadays has achieved the same status as a true observable Q. 

The increased accuracy of the data as well as the possibility of performing more elaborate nu- 
merical calculations, has stressed the need for better descriptions of (e,e'p) data. The source of the 
spectroscopic factor has been investigated following various hypotheses, from relativistic effects to 
in medium modifications of the nucleon properties [^. At the moment, one of the most frequently 
considered hypotheses is the partial occupation of the single particle levels produced by the short- 
range correlations 0. Various estimates of these occupation probabilities have been done both in 
nuclear matter ^ and in finite nuclei |ll|, but their values do not seem to be compatible with 
the empirical spectroscopic factors. On the other hand, direct calculations of (e,e'p) cross sections 
with models taking into account short-range correlations show noticeable effects produced by these 
correlations |12]-[15|, in spite of the fact that also in these cases spectroscopic factors are still needed. 
These calculations of the semi-inclusive cross section have been done with simplified treatments of the 



correlations. In some calculations only the lowest order cluster terms have been considered [12, 13 



In other calculations |14] the overlap wave function, which in the mean field model corresponds to the 
hole single particle wave function, has been extracted from correlated one-body density distributions 
evaluated via more or less sophisticated ground state calculations. 

There are two weak points in these calculations. A first one concerns the inconsistency of the 
various inputs of the calculations. Single particle wave functions and correlations are not linked by a 
unique nuclear hamiltonian as they should be. The second weak point is related to the fact that the 
proper normalization of the many-body wave function is not preserved. 

We have developed a model which takes into account all the diagrams containing a single correlation 
line [|l^]. Since the cluster expansion conserves the sum rules order by order [17|, the normalization 
of the wave function is maintained. This is immediately seen in the application of the model to the 
description of the ground state densities p^, 19 1. The validity of our model has been tested in nuclear 



matter by comparing our results with those obtained in a calculation that considers all the cluster 
terms [20|. The agreement between the two results is excellent, showing that, for relatively simple 
operators, like the charge operator, the first order approximation is reliable. 

The inputs of our model, single particle wave functions and correlations, are taken from Fermi 
Hypernetted Chain (FHNC) calculations of the ground state of closed shell nuclei obtained by mini- 
mizing the hamiltonian expectation value evaluated with realistic and semi-realistic nucleon-nucleon 
potentials |^, ^ ]. 

In this paper we present the results of our investigation of (e,e'p) reactions in ^^O. We have 
analyzed the validity of various approximations commonly adopted in the mean field description of 
these processes. It is not at all evident that approximations which are within a certain theoretical 
framework should remain valid in a different one. Moreover, if the aim is to disentangle subtle physical 
effects, the uncertainty in the result due to some adopted assumption can be comparable with the 
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magnitude of the searched effect. 

Summarizing in few words our findings we may say that short-range correlation effects are ex- 
tremely small and cannot explain the values of spectroscopic factors needed to reproduce the data. 

2 The cross section 

In this section we briefly recall the expressions of the coincidence cross section used in our calculations. 
We work in natural units {H = c = 1, e'^ = 1/137.04) and employ the conventions of Bjorken and Drell 
[p^]. The initial and final electron four vectors are respectively k = (e,k) and k' = (e',k'), and we 
used the symbol q = {uj, q) = k — k' to indicate the four momentum transfer. The four momentum of 
the emitted nucleon is p = (ep,p). The reference system, shown in Fig. 1, has been defined following 
the commonly used prescriptions. The scattering plane is defined by the electron vectors k and k', 9 
is the angle of the scattered electron, the quantization axis is taken along the direction of q, and the 
angle 9p and (pp define the direction of the emitted nucleon. 




Figure 1: Reference system used in our calculations. 

In our calculations the cross section has been derived by adopting the commonly used assumptions 
[p^]: the electron wave functions are plane wave solutions of the Dirac equation, only one photon is 
exchanged between the electron and the nucleus, all the terms depending upon the electron rest mass 
are neglected. With these approximations we obtain: 
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where we have indicated with cta/ the Mott cross section: 
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Since we describe the emitted particle within the non relativistic kinematics, we obtain K = |p|mp. 
The factors v come from the leptonic tensor and depend only from kinematic variables: 
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Vtt = -4- (6) 
2q2 ^ ^ 



The information about the nuclear structure is included in the w factors. Because of the current 
conservation only three components of the current four-vector are independent. We choose the charge 
p(q) and the two transverse components in spherical coordinates: 
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J±=T—{Jx±iJy). (7) 



We can express the w factors as [25 



Wl = {^,\pHclmf){^Mq)\^^), (8) 

Wt = (^i|jl(q)|vI/y)(%|J-(q)|^.) + (^i|4(q)|*/)(M//|J+(q)|^.), (9) 
Wtl = 2i?e((*i|^t(q)|^^)(xj/^|j_(q)|vl;.)_(^.|pt(q)|^^)(^.|jt(q)|^^^^ ^ (10) 

Wtt = 2iie((^,|j|(q)|^/)(^y|J_(q)|M/,)) , (11) 

where we have indicated with and |^/) initial and final states of the full hadronic system. In the 
previous equations a sum on a and mh the projection on the z axis of the spin of the emitted particle 
and of the angular momentum of the the residual nucleus it is understood. 

In this work we consider only one-body electromagnetic currents. The charge operator is expressed 

as: 

^ 3 

= E^^'^Cr-r.), (12) 

the convection current operator as: 

^ ■ 3 

•^""(r) = E ^ ['^(r - r.)V, - V,<5(r - r,)] , (13) 

and the magnetization current operator as: 

= t i ("^H^ + ""H^) V X --0 . (14) 

In the previous equations indicates the rest mass of k-th nucleon, /i^ and fi^ the anomalous 
magnetic moment of the proton and the neutron respectively, cr^ the Pauli spin matrix of the k-th 
nucleon and r| = 1 (—1) for protons (neutrons). In our calculations the nucleonic internal structure 
has been considered by folding the point-like responses with the electromagnetic nucleon form factors 
ofRef. 11. 



3 The nuclear model 

The nuclear final state \^f) is described in the asymptotic region as the product of the wave functions of 
the emitted nucleon, (j){p-r)xcT, with p its momentum, and of the rest nucleus, \I'(ri...rA-i; en] jh,iT^h)- 
We have indicated with the angular momentum of the rest nucleus and we have defined = 
ep — uj. In a strict independent particle model these quantities correspond to the energy and angular 
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momentum of the hole state. We perform a multipole expansion of \^f) to express it in terms of 
eigenstates of the total angular momentum of the A nucleon system: 

IpfJ-pjpnip JM,n 

{jp'mpjhmh\JM)\^!; JM, U; {Ipjpnipep, Ihjhmheh)) 
= SE E i'''yiUpiP)i^pf^pl'^\jp^p) 

IPI P JM,U ^ 

{jpmpjhmh\JM) l-^^' ^.p, h) 

{JM,U;p,h\JM,U;p, h)2 

In the above equations the state |^'; JAf, H; {Ipjpnipep, Ihjhmheh)) represents the excited state of 
the system with A nucleons with total angular momentum J projection M, parity 11, having a particle 
in the continuum wave characterized by orbital and total angular momentum Ip and jp with projection 
nip and energy ep, and a residual nucleus with hole quantum numbers Ih jh, and eh- We have 
indicated with Yj^ the spherical harmonics and with the symbol {lamalbmi)\J M) the Clebsch-Gordan 
coefficients [^]. In the last expression, the nuclear excited state is given in terms of unnormalized 
many-body wave function | JM, 11;^, h) 

In our model, the nuclear final state with total angular momentum J is described as a superposition 
of all possible one-particle one- hole transitions allowed by the angular momentum selection rules. 
Within this approximation, the knowledge of the energy and of the momentum of the emitted nucleon 
determines the quantum numbers of the hole state of the residual nucleus. Then the summation on 
the hole states drops and we sum only over the possible particle states in the continuum. We defined 
the sum over p as a summation over the Ip, fip, jp and rup quantum numbers. 

From the equations (P)-([Tl|) defining the w factors, one can see that the basic ingredients to be 
calculated are the transition matrix elements from the ground state to the excited state of the nuclear 
system induced by a one body operator. We write the transition matrix elements as: 



{^f\Or,{ci)\^i) = nE E i-iy'Ylpf^piP){^pf^pl^\jp^p){jp'n^pjhrnh\JM) 

(JM, n;p, /i| 0^(q) 1 00, +1) 



( JM, n; p, h\JM, n; p, h) 2 (00, +1|00, +1) 2 
An 

= ^E E "^(P'^S'^'^'n) ej,Af,„,n(q;p,/i), (16) 

IPI P J,M,^ 

where |00, +1) is the unnormalized ground state which, for the nuclei we shall consider, has zero 
angular momentum and positive parity. We have included all the geometric part of the expression in 
the symbol 

A{p,h;J,M,U) = {-iy''Yi^f,^{p){lpfip-a\jpmp){jpmpjhmh\JM) , (17) 

and we have indicated with 0,j(q) the generic one-body operator inducing the transition. For rj = 
this is the charge operator, while for rj = ±1 is given by the appropriate sum of convection and 
magnetization currents, Eq. (|^). 

In our nuclear model we consider the nuclear states described as: 

|00,+1) = F|$;00,+1), (18) 
\JM,U;p,h) = F \^;JM,U;p,h) . (19) 

We have indicated with |$; 00, +1) the Slater determinant describing the mean field wave function in a 
pure independent particle model. This means that, taken a basis of single particle wave functions, all 
the states below the Fermi surface are fully occupied and those above all completely empty. The state 
1$; JM, Il;ph) indicates a Slater determinant where the hole function h has been substituted with the 
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continuum particle function p. With respect to a pure independent particle model, the novelty is the 
presence of the correlation function F. This function has, in principle, a very complicated operatorial 
dependence, analogous to that of the hamiltonian. Our calculations have been restricted to the 
use of purely scalar (Jastrow) correlations. For this reason we immediately simplify the expressions 
formulating them only in terms of this type of correlations. The adopted ansatz on the correlation is: 

A 

Fil,2,...A) = l[fiu,), (20) 

i<j 

where rij = |rj — rj| is the distance between the positions of the particles i and j. 

Using the well-known cluster expansion techniques we rewrite the transition matrix 0,M,»?,n(q;Pi h) 
in Eq. (^) as: 



(JM, n;p,/i|O^(q)|00,+l) 
(00,+l|00,+l) 



(00,+l|00,+l) 



{JM,U;p,h\JM,U;p, h) 



(21) 



The two factors in Eq. ( [2l| ) are separately evaluated by expanding both numerator and denominator in 
powers of the short-range correlation function. The presence of the denominators is used to eliminate 
the unlinked diagrams |2£]. 

Since in our calculations the correlation functions are purely scalar they commute with the operator 
0^(q) and therefore we have to deal with terms of the kind: 

(JM, n;p,/i|O^(q)|00,+l) = JM, n;p, /i|FtO^(q)F|$; 00, 

A 

= JM, n;p, /i|0^(q) J] f\nj)\^- 00, 

i<j 
A 

= JM,n;p,/i|O^(q)n(l + /i*j)l^;00,+l)L, (22) 

where we have used the function hij = f'^{rij) — 1 and the subindex L indicates that only the linked 
diagrams are considered. 

The approximation of our model consists in retaining only those terms where the hij function 
appears only once: 

0,M,r7,n(q;p,/i) — > Cj,M,r,,n(q;p,/i) 

= ($;JM,n;p,/i|0^(q) ^(l + /i,,)|$;00,+l)L. (23) 



This result has been obtained using a procedure analogous to that adopted in Ref. [18| for the 
evaluation of the density distribution, and therefore the truncation of the expansion is done only after 
the elimination of the unlinked diagrams. 

We show in Fig. 2 the Mayer like diagrams describing all the terms considered in our calculations. 
This set of diagrams has already been presented elsewhere [16, but we show them again because 
the identification of the individual diagrams is essential for the discussion of the results. In each 
diagram the black square indicates the coordinate where the one-body operator 0^(q) is acting, while 
the black dots indicate the other coordinates. The dashed line represents the correlation function /ijj, 
which operates on two-coordinates only, and the continuous oriented lines the single particle wave 
functions 4>k- 

If we label with 1 the coordinate where the external operator 0(q) is acting, we can specify the 
of Eq. (ID as: 

eWn(q;P,/i) = JM, n; p,/i|0^(q) I 00, +1) 
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Figure 2: Mayer-like diagrams considered. The black squares represent the point where the external field is 
acting on, the dashed line the correlation function and the oriented lines the single particle wave functions. 



+ JM,n;p,/i|0^(q) ^ /i^ |$; 00, +1)^ 

A 

+ {^;JM,U;p,h\0^{q) ^ hij \^;00, +1) l ■ 

i<i<j 



(24) 



The above expression shows that our model generates, in addition to the uncorrelated transitions 
represented in Fig. 2 by the one-point diagram (1.1), also two- and three-point diagrams. The presence 
of these last diagrams is necessary to have the correct normalization of the many-body wave function, 



as discussed in 1161. 



In our case the operators 0,j(q) are the Fourier transforms of the charge and current operators 
0^(r) given in Eqs. (|T2|)-(p!^). Since we describe the nuclear transition between states with good 
angular momentum, it is convenient to make a multipole expansions of these operators. For the 
charge we have: 

MjMiq) = f d\jjiqr)YjM{r)p{r), (25) 



where we have indicated with r the and (f) angles characterizing the vector r in polar coordinates 
and with jj the spherical Bessel functions. 

For the current we should distinguish between the electric excitations (E), with natural parity 
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n = ( — 1)'^, and the magnetic excitations (M), with unnatural parity 11 = (—1)'^"'"^: 

1 



dV |v X 



JJ 



(gr)Y)S(?)]}-J(r) 



and 



T^M = |dVj,(gr)YjS(f-)-J(r), 



(26) 
(27) 



where we have used the symbol Yjj to indicate the vector spherical harmonics [p^ ]. 

By substituting the above equations into Eq. (jl^) and using the Wigner-Eckart theorem and the 
properties of the 3j symbols, we can rewrite Eq. (|16|) in terms of reduced matrix elements: 



(^^|0,(q)|^,) 



An 



J2 E AP^h;J,M,U) C\M,rj,ni^-,P,h) 



P J,M,U 



p J,n 



xi-iY 



(CD; J,n;p,/i||Oj(^)(q)||cI>;0,+l) 



+ ($; J,n;p,/i||Oj(^)(q) E \\<^;0, +1) 

i>i 
A 

+ (<!>; J, n;p, /ill Oj(^)(q) E \\<I>;0,+1) 

i<i<j 



(28) 



where we have defined the operators: 



(29) 



if ?7 = 
if ?7 = lb 

In Eq. (^) the sum on M drops because M = —rj due to angular momentum conservation. 

Following Ref. [16| the calculations of the transition matrix elements are carried out by performing 
a multipole expansion of the correlation function hij: 



hij = h{rij) = h{r i, rj, cos 9 ij) = E hLiri,rj) PL{cos6ij) 



(30) 



where Pl are the associated Legendre polynomials. 

The final state quantum numbers are determined by the quantum numbers of the uncorrelated 
many body state I^ja/) built as a Slater determinant of single particle wave functions of the form: 



1 



(31) 



where Xs and xt are the spin and isospin wave functions respectively. 

By using the orthonormality of the single particle basis, we rewrite Eq. (p^) as: 



{^f\Or,{cD\^^) 



47rV27r 



IPI 



p J,n 



Oj(^)iq)\\h) + J2iP^j\\Ojin)i^)hi2\\h,j) 



+ E (P'i'^ll Oj(^)(q) h23\\h,j,k) 

jk 

E E(-^)"^-^(P' J, -r?, n) ej,,,u{^; P, h) , 
IPI ph j,n 



(32) 



7 



where (,j^Yi{c[',p,h) is the sum of the reduced matrix elements corresponding to the one-, two- and 
three-points diagrams. Here the sums of the j and k indices are running on all the hole single particle 
wave functions. The subindices 1,2,3 of the correlation function h are associated with the radial 
coordinates of the first wave function, where also Oj(^)(q) is acting, the second and the third one 
respectively. We have indicated only the direct matrix elements in the second and third term of the 
equation, but we calculate all the exchange diagrams, as shown in Fig. 2. 

The multipoles of the one-body operators are calculated by inserting the expressions ([l^ ) 
and (H) into Eqs. (||), (|2|) and (|2|). Generally, one may write the resulting operator Ojj^j{q) as 
the product of a term depending on q and on the modulus of ri times a term depending only on the 
angular coordinates: 

0'jj,,{ci)^Fj{qn)OjM{ri). (33) 
The evaluation of the matrix elements of the various diagrams for the three operators follows 



the lines of Ref. |16]. In the present calculation we have added all the correlated diagrams of the 
convection current, which in that reference was included only at the mean field level. 

All the ingredients needed to calculate the cross section are now available. The information on the 
angular distribution of the emitted particle is fully contained in the spherical harmonics of Eq. (|l7|), 
which has to be squared to calculate the response. We have to deal with a geometrical term of the 
type: 



(^i-ir''Aip',h';J',-v',U')J (^( 



\-J 



A{p, h; J, -7], n) 



{Ipi-nih - ?/ - a)-a\jp{-mh - r]')){jp{-mh - r]')jhmh\J' - r]') 



{lp{-mh - 77 - (^)7;(^\jpi-mh - r])){jpi-mh - v)jhmh\J - v) 



{-it 



I'D Ir, J ~\~J 



]p,{-mh-ri-a) (p)^|],(-mh-»7'-a) (p) 



(34) 



This expression arises because the angular momentum of the hole jh , its z axis projection ruh, and a 
the third component of the spin of the emitted particle, are good quantum numbers of the hadronic 
final state. The above equation shows that the dependence from the angle (j)p between the scattering 
plane and the plane where the momentum of the emitted particle is lying, is present only when rj ^ r]', 
i.e. only in the interference terms tl and tt in Eq. (||). 
We rewrite the cross section expression as: 



a 



de'dVL.dE'd^r. 



(27r)3 



-2vuRe{Wo,-i - VFo,+i) cos (t)p + 2vttRe{W+i,^i) cos 2(Pp] 



(35) 



where the W functions are defined as 



327r3 



{lp{-mh -r]- a)^a\jp{-mh - 



,2 ■^h-r])){jp{-^h-v)jhmh\J -r]) 

{Ipi-nih -ri - a)^a\jp{-mh - 'n')){j'p{-mh - v')jhmh\J' - v') 

{lp + mh + rj + a)\ {I'p + nih + y' 
(1 - ^, - ^ - (V - ^, - ^1 



{Ip - rrih - rj - a)\ {I'p - rrih - rj' - a)! 
p-™''-''-"(cos Op) P^7"'^"'''""(cos Op) 

6,7?,n(q;p,/i) Cjt^',n'(q;f''' 



2lp + 1 



4tt 



21' +1 



,h). 



(36) 
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4 Specific applications 



We have applied our model to the investigation of the ^^0(e,e'p)^^N reaction in the quasi-elastic 
region. This nucleus has been chosen for two main reasons. The first one because it is a sufficiently 
light nucleus, such that the plane wave description of the electron wave functions is still valid [pC|]. 
The second reason is related to the fact that this nucleus has been studied with microscopic theories 
that provide us the needed input for our calculations. Specifically, the single particle bases and the 
correlation functions are taken from Ref. |21| and [22| where they have been fixed by minimizing 
microscopic many body hamiltonians. The calculations of Ref. ||2^ have been done with the semi- 
realistic nuclear interaction S3 of Afnan and Tang [^] by using a purely scalar correlation function. 
The calculations of Ref. [p^ j have been done with a reduction of the Argonne V18 interaction [^ ] 
called V8' and used to perform Green function Monte Carlo calculations |33|. In |22] a complicated 
state dependent correlation functions was used, while in the calculations of the present paper we have 
considered only the scalar term of this correlation. These two inputs are the same used in Ref. [^] 
to calculate the inclusive responses, and as in that case we shall label with S3 and V8 the results 
obtained respectively with them. For sake of completeness we give in Table |l] the parameters of the 
Woods-Saxon potentials and we show in Fig. 3 the two correlation functions used. If not specifically 
mentioned our results have been obtained by using the same mean field for particle and hole states. 
We shall investigate in Sect. the effect of using different mean fields as it is commonly done in the 
investigation of (e,e'p) data. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



r [fm] 

Figure 3: Correlation functions used in our calculations. 



Table 1: Parameters of the Woods-Saxon potentials used to generate the set of single particle wave 
functions. 





V8 S3 


vr 


Vq [MeV] 
Vis [MeV] 
R [fm] 
a [fm] 


-53.0 -52.5 
0.0 -7.0 
3.45 3.2 
0.7 0.53 




Vo [MeV] 
Vis [MeV] 
R [fm] 
a [fm] 


-53.0 -52.5 
0.0 -6.54 
3.45 3.2 
0.7 0.53 



Our calculations have been done in both perpendicular and parallel kinematics. In the first sit- 
uation the value of the momentum transfer q is kept fixed and the changes oi pi = |p — q|, the 
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momentum of the proton before being emitted, are obtained by changing the detection angle 9p. This 
is, for us, the easiest case to calculate and to investigate, even though all the four response functions 
Wj^rji contribute. In the second situation, the parallel kinematics, the momentum transfer q and the 
momentum of the emitted nucleon p are kept in the same direction and the missing momentum is 
changed by changing the values of q. In this case the interference responses, those with t] ^ r]' do not 
contribute p5[| . 

Since in (e,e'p) processes the number of variables into play is relatively large, we have mainly done 
our investigation at a fixed value of to = €i — ej = 128 MeV, and for a fixed value of momentum of 
the emitted proton, |p|= 444 MeV/c, knocked out from the lpi/2 level. In perpendicular kinematics 
the scattering angle was fixed to set q = p, to be able to scan all the values of pi from up to 2q. 
The kinematical conditions just described correspond, for the perpendicular kinematics, to those of 



experimental data of Ref. |34| and, for the parallel one, to the data of [p5| 





Figure 4: Cross sections of the "'^^0(e,e'p)^^N process for the kinematics described in the text as a function of 
the nucleon initial momentum pi = \pi\ = |p — q|- The proton is emitted from the lpi/2 level. The left panels 
have been obtained in perpendicular kinematics while the other ones correspond to parallel kinematics. The 
labels S3 and V8 refer to the two different inputs as discussed in the text. The full lines represent the results 
obtained with our model, while the dashed ones have been obtained leaving out the three-body diagrams of Fig. 
2. The mean-field results are almost exactly overlapped to the full ones and are not included in the figure. 



4.1 General features of the results 

In this section we discuss the general features of our results. In Fig. 4 we present the cross sections for 
the kinematics described above as a function oipi. We should recall that we defined = \pi\ = |p — q|. 
In perpendicular kinematics, we consider the values of pi to be positive when (pp=0 and negative when 
0^=180*^. In parallel kinematics pi > when |p| > |q| and pi < when |p| < |q|. In Fig. 4, the 
left panels show the perpendicular kinematics results while the right ones correspond to the parallel 
kinematics. The full lines have been obtained considering the complete model. Since the main novelty 
of our calculations is the presence of the three-body diagrams, we shall discuss their role in some 
detail. For this reason we show with dashed lines the results obtained by adding to the mean field 
only those terms coming from the two-point diagrams. The mean-field results we have obtained, are 
almost exactly overlapped to the full ones, and have not been shown. 
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The first remark we can make is that the effect of the correlations is very small, as it was found in 
Refs. 36] for the inclusive reactions. As already pointed out in Refs. [0, |l^, |l^ the three-point 
diagrams reduce the effect produced by the two-point diagrams whatever it is. This is the consequence 
of the requirement of the conservation of the number of particles. 
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Figure 5: Normalized differences between the correlated and uncorrelated cross sections of Fig. 4. These are 
the differences between the cross sections divided by the maximum values of the uncorrelated cross sections. 
The dashed lines represent the two-point diagrams results while the full lines those obtained with the full 
calculations. The left panels show the perpendicular kinematics results and the right panels the parallel ones. 



In order to emphasize the effect of the correlations we present in Fig. 5 the same results in terms of 
difference between the correlated cross sections and the uncorrelated ones. In order to have a relative 
measure of the effect we have divided these differences by the maximum values of the uncorrelated 
cross sections. Henceforth we shall call these results normalized differences. 

The effect of the two-point diagrams is rather different for the S3 and V8 calculations. In the 
former case is positive, the mean field cross section is increased by the correlations, while it is negative 
in the latter case. It is interesting to notice that the final result is rather similar in both cases as one 
would expect given the strong similarity of the two correlation functions used. 

The effect of the correlations on the individual responses calculated in perpendicular kinematics 
is shown in Fig. 6. We should first notice that the charge response Wl = VFo,o and the transverse 
response Wt = + W-i-i are much larger than the interference responses Wtl = 2i?e(T^_i,o — 

VFo,-i-i) and Wtt = 2iZe(VF+i^__i). In the charge response Wl the various terms of the correlations 
behave like the full cross sections, the final results being a lowering of the mean-field response. This is 
not the case for the Wt response, showing an increase of the mean-field results. This effect has been 
discussed in detail in Ref.[|16| where it has been shown that it is related to the absence of the diagram 
(2.3) of Fig. 2 in the longitudinal response. 

In what refers to the interference responses, the behavior of Wtl is similar to that of Wl- In Wtt 
the correlations increases the mean field responses. In this last responses the 3-point diagrams appear 
to be very sensitive to the input, their contribution is negligible for V8 calculation and noticeable for 
the S3 ones. 
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Figure 6: Differences between correlated and uncorrelated response functions in the perpendicular kinematics 
for the calculations of Fig. 4. The differences have been divided by the maximum value of the uncorrelated 
responses. The full lines represent the results of the complete calculations and the dashed ones those obtained 
using the two-point diagrams only. 
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Figure 7: Normalized differences between correlated and uncorrelated cross sections calculated in perpendicular 

kinematics for different values of the momontuni transfer q = |q|. As in the previous figures the dashed lines 
refer to the two-point diagrams calculations while the full lines show the results of the full calculation. 
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The cross section is obtained by summing the contributions of all the responses. Prom Eq. (35) one 
can observe that the Wtl responses is added for positive values of pi and subtracted for the negative 
ones, and this generates the asymmetry in the perpendicular cross sections of Fig. 4. 

We have investigated the momentum transfer dependence of the correlation effects by calculating 
(e,e'p) cross sections in perpendicular kinematics for different values of |q|. In all the calculations 
we have fixed the momentum of the emitted nucleon such that |p| = |q|. Results obtained for four 
different values of |q| are shown in Fig. 7 as normalized differences. As in the previous figures the 
dashed lines show the results obtained by adding to the mean field calculation the two-point diagrams 
only, while the full lines show the results of the full calculation where also the three-point diagrams 
have been added. The shape of the lines is slightly modified at different values of |q|, but the order 
of magnitude of the effect is almost the same. It is interesting to notice that the cancellation of the 
two-point diagrams effects produced by the inclusion of the three-point terms generates results rather 
independent from the correlation function. 
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Figure 8: Cross sections for the emission of the IP3/2 proton in perpendicular kinematics for the S3 and V8 input 
parameters. The dashed Hnes show the results obtained adding the two-point diagrams to the uncorrelated terms 
and the full lines the result of the complete calculation. The uncorrelated results are almost exactly overlapped 
to the full ones and have not been shown. In the panels a and b we show the cross sections while in the panels 
c and d the normalized differences between correlated and uncorrelated cross sections. 



We have also verified the possible changes of the correlation effects depending upon the hole wave 
function. As an example we show in Figs. 8 and 9 the cross sections, and the normalized differences, 
obtained for the emission of the IP3/2 and lsi/2 protons under the same perpendicular kinematics 
of Fig. 4. While for the IP3/2 case the shapes of the cross sections are noticeably similar to those 
previously shown, the curves obtained for the lsi/2 proton are rather different. The small differences 
in the radial shapes of the two p hole wave functions do not modify sensitively the behavior of the cross 
section. On the other hand, the big difference in the s wave cross section is due to the different radial 
shape of the hole wave function. In spite of this the behavior of the correlation effects is analogous 
in all the cases we have considered. The calculations done with two-point diagrams show noticeable 
sensitivity to the change of the correlation, while, when the three-point diagrams are added all the 
results are similar and almost independent of the correlation function. 
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Figure 9: The same as the previous figure for the emission of the lsi/2 proton 



4.2 Orthogonality and consistency of the input 

In the analysis of the (e,e'p) experiments it is a common practice to use a real mean field potential 
to describe the hole wave function while the particle wave function is described using an optical 
potential containing also an imaginary part. This procedure, adopted to take into account the final 
state interaction of the emitted particle with the rest nucleus, produces a non orthogonality between 
particle and hole single particle wave functions. The effects of this procedure on the (e,e'p) cross 



section have been thoroughly studied in Ref. |37] where it has been shown that they are negligible for 
the typical kinematics used in the experiments. 

We have investigated the non orthogonality effects when correlations are included in the calcula- 
tions. In our derivation of the transition matrix elements the orthogonality of the single particle wave 
functions is essential. On the other hand, since we use central potentials, only the radial parts of the 
wave functions change, therefore all the angular momentum algebra is still valid. This means that 
the non orthogonality is present only between single particle wave functions having the same angular 
momentum but different number of nodes in the radial part. 

We have done two calculations where the particle wave functions have been evaluated with a 
different mean field with respect to the hole wave functions. In the first one the particle mean field 
was set to zero; therefore the emitted particle wave function was described as a plane wave. This kind 
of calculation is commonly referred as the Plane Wave Impulse Approximation (PWIA). In the second 
type of calculation we used the optical potential of Ref. [^] with the parameterization adapted for 
in Ref. ||]. 

We show in Fig. 10 the cross sections calculated with the S3 and V8 inputs for the same per- 
pendicular kinematics conditions used in Fig. 4. The panels a and b show the PWIA results, while 
the other two panels those obtained with the optical potential. A comparison with Fig. 4 shows that 
the PWIA cross sections are larger than the original ones, while the cross sections obtained with the 
optical potential are smaller. It is interesting to notice that the PWIA results have a minimum for 
Pi=0 which is exactly zero, at expected, because in this approximation the cross section is related to 
the Fourier transform of the hole wave function. 
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Figure 10: Cross sections calculated in perpendicular kinematics under the same conditions of Fig. 4 with 
different mean field potentials for the particle wave functions. In panels a and b the results obtained in Plane 
Wave Impulse Approximations are presented. In panels c and d those obtained the optical potential of Schwandt 
et al. pa. The meaning of the line is the same as in Fig. 4. 



In both kinds of calculations the effect of the correlations is very small and to emphasize it we 
show in Fig. 11 the normalized differences. The behavior of the various terms is similar to that shown 
in Fig. 5. The results obtained with the two-point terms are rather different when calculated with 
S3 and V8 inputs. When also the three-point diagrams are included they become similar. Clearly 
the effects of the non-orthogonality of the basis are negligible on the correlations. We have obtained 
similar results also in parallel kinematics. 

A comparison between the three calculations done in perpendicular kinematics shows that the 
shapes of the results obtained in PWIA and with the optical potential are rather similar and quite 
different from the results obtained with the same real potential used for the hole states. This arises 
because the real part of the Scwhandt potential is quite weak, —28.359 MeV, against the about —50 
MeV of the potential used in our V8 and S3 calculations. The effect of the optical potential is to 
conserve the PWIA shapes but to lower the values of the cross section. This second effect is obtained 
because of the presence of the imaginary part of the potential. 

We have already stressed that our inputs are taken from microscopic FHNC calculations where, 
for a given hamiltonian, the energy mean value has been minimized by varying both mean field and 
correlation function. This consistency of the input is not usually respected in (e,e'p) calculations with 
short-range correlations, where mean field and correlation functions are taken from different sources. 

We have tested the need of relating mean-field and correlation functions by interchanging them in 
our inputs. This means that we have done calculations using the S3 mean field potential together with 
the V8 correlation functions, and vice versa. The results of these calculations are shown in Fig. 12 as 
normalized differences and are compared in each panel with the results of Fig. 5. Also in this figure 
the need of including three-point diagrams for the reasons already mentioned becomes evident. The 
final results seems to be rather independent of the mean field used. On the other hand, the quantity 
shown has been defined to minimize the effects of the mean field and enhance those of the correlations. 
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Figure 11: Normalized differences between the cross sections presented in Fig. 10. The meaning of the lines is 
the same as in Fig. 5. 
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Figure 12: Normalized differences between cross sections calculated with the same mean field potential but 
with different correlation functions. In panels a and b the S3 mean field has been used, while in panels c and d 
the V8 mean field was considered. The S3 and V8 labels indicate the correlation function used. The calculations 
of the a and c panels have been done in perpendicular kinematics and those of the b and d panels in parallel 
kinematics. As usual the dashed lines refer to calculations with two-point diagrams while the full lines represent 
the complete calculations. 
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4.3 Comparison with other approaches 



In this section we shall try to make a comparison between our approach and other models that handle 
short-range correlations in electron scattering. Depending upon the methodology used to attack the 



problem we classify these models into two categories. A first kind of calculations |15] uses an approach 
similar to ours but the evaluation of the diagrams remains restricted to the so-called single pair 
approximation. This approximation implies that only those two-point diagrams where the knocked 
out particle is directly connected to the electromagnetic operators are considered. In our classification 
scheme these are the diagrams (2.1) and (2.2) of Fig. 2, even though the exchange diagram (2.2) is 
not always considered. 

The second category of calculations ||l^, Q is based upon the idea of the spectral function |Q, ^. 
The basic hypothesis is that the nuclear finite state wave function can be separated in terms of emitted 
particle wave function and spectral function which is related to the overlap between the wave function 
of the target nucleus and that of the A-1 nucleon system. Spectral functions used in the evaluation 
of (e,e'p), and also (e,e'2N), cross sections have been calculated in various manners. In some case the 



Brueckner G- matrix theory has been used [10|, in others they have been obtained from variational 
Monte Carlo calculations [^ ]. In the work of Ref. fl^ the spectral functions have been obtained 
using the asymptotic properties of the one-body density matrices obtained from microscopic theories. 
A precise comparison with our approach is rather difficult since the spectral functions include high 
order correlation terms which we do not consider. What remains, however, is the fact that also 
this approach considers only those diagrams where the emitted particle is directly connected to the 
electromagnetic operators. In our approach this could be translated by saying that in addition to the 
diagrams (2.1) and (2.2), already mentioned, also the diagrams (3.3) and (3.4) are considered. 

We have tried to mimic the other approaches by repeating the calculations of Fig. 4 without 
those diagrams where the particle is not directly connected to the electromagnetic operators. In these 
calculations the 1/2 and 1/3 terms multiplying the two and three point diagrams respectively have 
been set to 1, as is done in the literature. 
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Figure 13: Normalized differences between correlated and uncorrelated cross sections. The full lines represents 
the complete calculations, the dotted lines have been evaluated by adding to the mean field the contribution of 
the diagram (2.1) of Fig. 2, the dashed dotted lines, the contribution of the diagrams (2.1) and (2.2) and the 
dashed hnes, the diagrams (2.1), (2.2), (3.4) and (3.5). 
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The results of these calculations are presented in Fig. 13 as normalized differences. The full lines 
show the results of the complete calculations, those of Fig. 4. The dotted lines have been obtained by 
using the diagram (2.1) only, the dashed dotted lines have been obtained by adding also the diagram 
(2.2) and the dashed lines by adding the diagrams (3.3) and (3.4). 

The inclusion of the diagrams (3.3) and (3.4) changes radically the effect of the correlations, 
modifying shapes and signs. The final result in the case of S3 is rather similar to that of the complete 
calculation, but for the V8 input there is a noticeable difference. This instability of these results 
shows, in a clear manner, the need of including all the terms of the expansion. 



4.4 Comparison with experimental data 



In this section we compare the results of our calculations with some experimental data available on 
^^O. In the literature these data are presented in terms of a reduced cross section obtained by dividing 
the cross section by the electron- nucleon cross section and the kinematic factor K of Eq. (|^). In our 
calculations we have used the electron- nucleon cjcci cross section of Ref . [|^] . 



In Fig. 14 we make the comparison with the lpi/2 proton emission data of Ref. |34] measured 
in perpendicular kinematics. The full lines represent the complete calculations of Fig. 4, where all 
the diagrams are considered and the same mean-field potential for both particles and hole has been 
used. In the panels a and b the results for the S3 and V8 inputs are shown respectively. In the panels 
c and d the same curves have been multiplied by reduction factors to improve the agreement with 
the experimental data. These factors are 0.58 for the S3 result and 0.5 for the V8 one. Even after 
quenching the curves the agreement with the data is rather poor. 
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Figure 14: Reduced cross section compared with the experimental data of Saclay |34 . The full lines have been 
obtained from the complete calculations of Fig. 4 in perpendicular kinematics. The dashed lines have been 
calculated with the optical potential of Ref. |Q. In panels c and d the curves of the upper panels have been 
multiplied by reduction factors to reproduce the data (see text). 



The dashed lines of the figure have been obtained by changing the particle mean field with the 
optical potential of Schwand et al. |38|, in the parameterization of Ref. [35|. Also in this case a 
quenching of the curves is necessary, and the final results, shown in the two lower panels, have been 
obtained multiplying by 0.58 and 0.5 the S3 and V8 results respectively. The dashed curves show a 
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better agreement with the experimental data. 
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Figure 15: Reduced cross sections compared with the experimental data of Nikhef |3^. The hill lines have 
been obtained from the complete calculations of Fig. 4 in parallel kinematics. The dashed lines have been 
calculated with the optical potential of Ref. p^ . In panels c and d the curves of the upper panels have been 
multiplied by a common reduction factors of 0.8. 

An analogous situation is found in Fig. 15 where the comparison is done with the data of Ref. 
[^ ], again for the emission of the lpi/2 proton but measured in parallel kinematics. The meaning 
of the curves is the same of the previous figure, and in this case we used for all the curves the same 
quenching factor of 0.8. 

The results presented in Figs. 14 and 15 clearly show that the fully consistent calculations are 
unable to reproduce the data. Not only in terms of their magnitude, but the obtained results are 
qualitative different with respect to the data distribution. The mean fields fixed by the minimization 
process in Ref. |21] and |22| are too strong to describe properly the motion of the emitted particle in 
the continuum as we have already observed. The shift of the minimum of the reduced response with 
respect to pj = in parallel kinematics is an effect produced by the distortion potential. 

In Fig. 16 we compare the results of our calculations with the empirical values of the electromag- 



netic responses given in Ref. |41|. These calculations have been done with the S3 input and using 
the Scwhandt et al. optical potential ||3^. The meaning of the various lines is that of Fig. 4: full 
lines show the complete calculations and dashed lines have been obtained by adding the two-points 
diagrams to the mean-field terms. We do not draw the mean-field results because they are almost 
overlapped with the full ones. All the lines have been multiplied by a reduction factor of 0.7. The effect 
of the short range correlations is very small. One can clearly identify the increasing of the responses 
produced by the inclusion of only the two-point diagrams, but the final results are essentially back 
to mean- field result. The calculation done with the V8 input produces analogous results, with the 



obvious difference in the two-point calculations. The same data have been studied in Ref. |42| where 
it has been shown that the effects of the Meson Exchange Current are larger than those we find for 
the short-range correlations. 
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Figure 16: Response functions compared with the experimental data of Ref. |Q. Our calculations have 
been done with the S3 input and using the optical potential of Ref. The full lines show the results of 

the complete calculations, the dashed ones those obtained with the two-point diagrams, and the dotted lines 
(overlapped to the full ones), the mean-filed results. The panel labelled Wl shows actually Wl + \ q2^^2 Wtt- 

The contribution of Wtt is in any case small. The response Wtl is multiplied by l/\/2 in order to agree with 
the response experimentally separated. All the curves have been multiplied by a 0.7 factor. 



5 Summary and conclusions 

In this work we have studied the effects of short-range correlations on the one-nucleon emission pro- 
cesses induced by electromagnetic probes. The description of these processes has been done by using a 
nuclear model that considers all the linked diagrams containing a single correlation line. This implies 
that in addition to the two-point diagrams, commonly considered in the literature, also three-point 
diagrams should be included in order to maintain the proper normalization of the many-body wave 
function. 

The calculations have been done using single particle wave functions and correlations taken from 
microscopic FHNC calculations We have considered only the scalar term of the correlation. 

We have calculated the ^^O (e,e'p) ^^N reaction for the perpendicular and parallel kinematics of the 



Saclay and Nikhef experiments of Refs. |34] and ||3| 

In the kinematic conditions investigated, the effects of the short-range correlations are extremely 
small as compared to the mean-field results. They are within the accuracy of the experimental data 
and the uncertainty in the theoretical inputs. 

The two correlation functions used in our calculations are rather similar, and therefore we expect 
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them to produce similar effects on the mean-field results. This does not happens if only two-point 
diagrams are included. The S3 increases the mean-field responses, and cross sections, while the V8 
lowers their values. The effects of the correlations become similar when the three-point diagrams are 
included. The detailed study of the separated responses shows that the three-point diagrams reduce 
the effect produced by the two-point diagrams alone, whatever they are. This fact is a consequence of 
the need of the three-point diagrams to ensure the proper normalization of the wave function. This 
result is almost independent of the momentum transfer or hole wave function. 

We have been concerned about the consistency of the input, and we have calculated cross sections 
with different mean fields describing the emitted particle, and we have also interchanged correlation 
functions and mean field in our inputs. We did not find significant changes on the effects of the 
correlations, saying that these effects are rather independent of the chosen mean field. 

We have tried to compare our short-range correlation model with other approaches one finds in the 
literature. In these approaches only those diagrams where the emitted particle is directly connected 
to the one-body electromagnetic operator are considered. We have done calculations where only these 
diagrams have been considered, and we found a strong dependence of the results on the number of 
diagrams included, showing big instabilities to small changes of the input. The situation is rather 
complicated by the fact that in these approaches the particle wave function is usually described as 
moving in an optical potential. Therefore it is not clear whether the diagrams we have calculated 
explicitly are effectively considered by the optical potential. 

The comparison with the experimental data of Refs. p4| , 35, 41 1, shows that the short-range 
correlation effects do not significantly improve the agreement with the data. In this respect vari- 
ous problems remain open, such as the physical origin of the required quenching or the reason why 
the parameters of the optical potential are so different from those describing the hole states. The 
investigation of these problems is beyond the scope of the present article. 

Finally we would like to discuss some of the limitations of our calculations related to the hypotheses 
and the approximations of our model. One of these limitations is due to the fact that we consider 
only those diagrams containing a single correlation function. Clearly higher order terms could be 
important. We have already mentioned that the test of the validity of our model was done in Ref. [^] 
for the inclusive nuclear matter charge response. The excellent agreement between our results and 
those obtained considering an infinite set of diagrams, give us confidence about our approximation, and 
we trust its validity also in the case of finite nuclear systems and for the one-body current operators. 
It would be interesting, in any case, to make other tests on it. 

Another limitation of our approach is related to the use of correlation functions of purely Jastrow 
type. Recent calculations show that the tensor correlations on two-nucleon emission processes |l^, ^] 
produce effect comparable, or even larger, than those induced by the scalar correlations. However, 
these calculations have been done in the single pair approximation, and we have shown how dangerous 
is to rely on such an approximation which can produce large effects since the three-point counter 
terms are not present. On the other hand, we know from microscopic calculations |44] that the tensor 
correlations are essential to bind nuclei. They are also important, together with Meson Exchange 
Currents, in the nuclear matter transverse quasi-elastic response and they are necessary to obtain 



a proper description of the electromagnetic responses of few-body systems |46]. The possible effect of 
the tensor correlations on nuclear responses is an open problem which should be investigated. 

Last, but not least, our model does not take into account those processes beyond the mean-field 
description of the nucleus falling under the generic name of long-range correlations. These correlations 
are related to the coupling between the single particle dynamics and the collective excitation modes of 
the nucleus. Their effects on (e,e'p) processes have been considered in a nuclear matter approach [pX|], 
but they do not seem to be relevant. However, one should consider that collective surface vibrations 
are not present in nuclear matter calculations, and they can be extremely important in finite nuclear 
systems, as an investigation on the charge density distribution differences of heavy isotones has recently 
pointed out |47]. 
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